%  演示调用 追赶法 分解
clear;
A=[2 -1 0 0; -1 2 -1 0; 0 -1 2 -1; 0 0 -1 2];   %对角占优的三对角矩阵
f=[1;0;0;0];

% 追赶分解法
x = my_forward_backward(A, f)

% 直接调用matlab左除运算符
x_ = A \ f